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Abstract 

The eigenvalue statistics of quantum ideal gases with single particle 
energies e„ = n a are studied. A recursion relation for the partition 
function allows to calculate the mean density of states from the asymp- 
totic expansion for the single particle density. For integer a > 1 one 
expects and finds number theoretic degeneracies and deviations from 
the Poissonian spacing distribution. By semiclassical arguments, the 
length spectrum of the classical system is shown to be related to sums 
of integers to the power a/(a — 1). In particular, for a = 3/2, the pe- 
riodic orbits are related to sums of cubes, for which one again expects 
number theoretic degeneracies, with consequences for the two point 
correlation function. 

1 Introduction 

Most investigations of quantum chaos have focussed on the effects in single 
particle systems. The prime examples of frequently studied systems, such as 
hydrogen in a magnetic field, the standard map or small molecules all belong 
to this class (Eckhardt 1988, Casati and Chirikov 1995). Even electrons in 
a solid, a standard many body system, has until recently been reduced to a 
single (quasi) particle system. Yet the study of many-body quantum systems, 
even if they are integrable can be interesting for several reasons. 

For once, the spin-statistics theorem, which requires quantum wave func- 
tions to be either totally symmetric (Bose-Einstein statistics) or totally anti- 
symmetric (Fermi-Dirac statistics) under exchange of particles changes the 
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spectrum compared to the simple Maxwell-Boltzmann type superposition 
of the individual single particle density. This gives rise to changes in the 
total density of states, as is well known in statistical mechanics. The full 
program of implementing the permutation symmetry semiclassically in the 
mean density of states and in the trace formulas has recently been taken up 
by Weidenmiiller et al. (1993a,b). 

Experimentally, small clusters are examples of systems with many de- 
grees of freedom for which one should use microcanonical averaging rather 
than canonical, since the number of atoms or electrons is rather fixed. The 
difference between microcanonical and canonical ensembles can be observed. 
The spectral statistics of small systems influences their thermodynamical 
behavior (Miihlschlegel 1991). 

Furthermore, symmetric systems can serve as a reference point for sys- 
tems with weakly broken symmetries. In particular, it has been observed 
that weak symmetries, perhaps of dynamical origin, can give rise to strong 
degeneracies in spectra, the Shnirelman peak (Shnirelman 1993, Chirikov 
and Shepelyansky 1995). 

In this contribution only integrable ideal gases will be analyzed. The 
energy levels are constrained to be a power of the quantum number, E n = 
n a , where n = 1,2,3,... and a > 1, except where noted. This family 
includes for instance the eigenenergies for a particle of mass m confined to a 
1-d box of width L, measured in units of 4tt 2 H 2 /2mL 2 . It also describes the 
asymptotic eigenvalues for particles in homogeneous potentials of degree p, 
for which a = 2p/{p + 1) (Seligman et al. 1985, Seligman and Verbaarschot 
1987). The harmonic case, a = 1, causes problems for the stationary phase 
approximations used below and will only be considered occasionally. The 
deviations of its level spacing distribution from Poissonian have been studied 
at length previously (Berry and Tabor 1977a). The quantum ideal gases then 
have the eigenvalues 



where the different statistics have been indicated. 

The questions addressed here concern the mean density of states (section 
2), the level spacing distribution in particular for integer a (section 3) and 
the behaviour of the pair correlation function (section 4). The results on 
the mean density of states are of a more general nature, whereas some ex- 
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amples are particular to the powers a considered. The final section contains 
some speculations on the relevance of number theory especially for the pair 
correlation function for some rational values of a. 

2 The partition function and the mean density of 
states 

2.1 The connection 

The partition function for a quantum mechanical system is defined as 

Z(/?)=tr e-P H = Y,e-P E i , (2) 

3 

where the last sum extends over all eigenvalues Ej. The density of states 
p(E) is related to Z by a Laplace transform, so that the poles in an asymp- 
totic expansion of Z for small (3 are related to the rate of divergence of p(E) 
for large energy and thus to te mean density of states. In particular, if 

z(/?)~5>/r 7i (3) 

3 

then 

3 ^ 

As an example, in the harmonic case a = 1, one has 

and for the mean density of states 

p(E) = l. (6) 

Thus, whereas the partition function contains contributions from positive 
powers of (3 (which contain information on the shortest periodic orbits in 
the system, Berry and Howls 1994), these terms are cancelled in the density 
of states by the poles of the T-function in the denominator. 
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2.2 The recursion relation for D-particle partition functions 

For the Maxwell-Boltzmann case with no restrictions on the integer sums, 
the partition function for D particles can be written down explicitly, 

Z^ B \0) = (Z 1 (/3)) D . (7) 

This is no longer possible for the symmetry reduced subspaces. However, 
there is a simple recursion relation involving the partition functions for all 
particle numbers up to D, 

z d(P) = ^J2(±i) k+1 z 1 (k(3)z D - k ((3) ; (8) 

fc=l 

the +1 applies in the Bose-Einstein subspace and the —1 in the Fermi-Dirac 
subspace. A combinatorial proof of this relation was given by Bormann and 
Franke (1993). A more direct analytical proof may be based on the grand 
canonical formalism (Reif 1965). 

The grand canonical partition function is defined as 

oo 

(l(z,/3) = J2 z ° Z d(P) (9) 

D=0 

oo 

= E E z D e-^n?, (10 ) 

D=Qn\,...,n D 

the last sum on quantum numbers is restricted by the selection rules for the 
different statistics. Passing to an occupation number representation, where 
g n denotes the number of particles in quantum state n, one obtains 

n(z,P)= ]T ^^"e-^n^". (11) 

91,92 — 

In the case of the Bose-Einstein statistics, the occupation numbers can take 
on all nonnegative integer values, and the summation on g n gives rise to 
geometric series. In the case of the Fermi-Dirac statistics, the occupation 
numbers can only take on the values and 1. The result for both cases can 
be combined in a single expression, 

OO _q 

n(z,P) = -H{l-zee-^ a ) , (12) 

n=l 
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where 9 = +1 for Bose-Einstein and 9 = — 1 for Fermi-Dirac statistics. 

The D-particle partition function follows from taking an ra-th derivative 
of 0, with respect to z at z = 0. The first derivative becomes 

71=1 

= S(z,/3)n(z,f3). (14) 

The recursion relation (||]) now follows from the formula for derivatives of a 
product and the observation that the k-th derivative of S at z = is related 
to the single particle partition function, 



d k S(z,(3) 



= 9 k Z x ((k + l)f3). (15) 

z=0 



dz k 

This result holds for other forms of the single particle energies as well. 
2.3 Examples 

The first few partition functions when reduced to the single particle partition 
function become 



MP) = \(Z 2 1 {0)±Z 1 



2 

Zm = ^(zf((3)±3Z 1 ((3)Z 1 (2(3) + 2Z l {3(5)) (16) 

Zm = ^ (Zf((3) ± 6Z 1 2 (/3)Z 1 (2/?) + 8Zx 09)^(3/3) 
+3Zf(2p) ±6Zi (4/3)' 

etc. 

where the plus signs apply for the BE-statistics and the minus signs for the 
FD-statistics. The corresponding densities of states become 

ME) = l/Z ni,n 2 5{E - n? - n?) ± i £ ~ 



p 3 (E) = I zZ 8{E-n?-r%-r%) (17) 



± I Y, S(E - n? - 2n?) + W S (E- 3n<*) 

ni,ri2 n 



etc. 
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The leading order contribution is Z^ MB ^ /D\, with the familiar permutation 
factor D\. Thus the mean density of states is to leading order the same 
for Bose-Einstein and Fermi-Dirac statistics and equal to 1/D\ the density 
of states for the Maxwell-Boltzmann statistics. The next to leading order 
terms correct for the counting of states with two or more identical quantum 
numbers. 

Further progress can be made if either the full single particle partition 
function or at least its asymptotic expansion are known. Such is the case 
for a = 1, the harmonic oscillator, and for a = 2 and n = 0, ±1, ±2, ±3, . . ., 
i.e. a particle on a ring. 

The single particle partition function for the harmonic energies a = 1 
and its asymptotic expansion have been given above (|5|). The partition 
functions for two particles then become 



*i E W = ^-i^ + Tiio^- (18) 

1 3 11 7 



For the mean density of states, these expansions imply 

P% E (P) = \E-\ (20) 
-P F 2 D {{3) = \E-\. (21) 

Here the correction terms have the same sign but different magnitude. 

As a second example we take a particle on a ring (periodic boundary 
conditions on the interval), for which 



zm= E e-/3 ~ i ■ ( 22 ) 



From this one can derive that 



Eigenvalue statistics in quantum ideal gases 



7 




and the densities 

UE) = \±^E-^ 

ME) = 1 E W±^* + ^E-W (24) 

Again, the plus sign refers to the Bose-Einstein statistics and the minus 
sign to the Fermi-Dirac statistics. Evidently, since the one particle partition 
function consists of the the leading order term only, all lower order terms 
are due to symmetrization and the signs clearly reflect the statistics. 

The size of the next to leading order corrections increase rapidly with the 
particle number D and become very important for large D. As can easily 
been shown by induction, 

D/2 

z D m = V r (D/2) i 1 ± CD(3l/2 ± • • •) (25) 

so that the mean density of states becomes 

^« = 5i^2) £(D - 2,/2 ( 1± ^- 1/2± -)- < 26 > 

with 

V2D(D - l)V{D/2) 1 5/2 

CD = V£T((z?-i)/2) ~W^ {D ~ l)/ ■ (27) 

The energy E c , where cr>E c ~ 1 thus increases like (D — l) 5 . This 
may be compared to the groundstate energy of the Fermi system, Ep ~ 
. Even when compared to this the importance of this term increases like 
E c /Ep ~ (D — l) 5 / 3 . The approach to the density of states of the classical 
ideal gas is thus very slow, and it will be difficult to estimate ground states 
accurately from the leading order terms (something that works surprisingly 
well in many cases in few degree of freedom systems). Some examples for 
the behaviour of the next to leading order corrections are shown in Fig. |]. 
Similar behavior can be found in the harmonic oscillator case. 
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Figure 1: Mean density of states for D-particle systems in the Bose (top) 
and Fermi (bottom) subspaces for particles on a ring. The densities are 
normalized by the desymmetrized Maxwell-Boltzmann densities. 

2.4 Asymptotic expansions of single particle partition func- 
tions 

For integer a one can use the Euler-MacLaurin summation formula to derive 
an asymptotic expansion for the single particle partition function starting 
from the representation 

^ 6 ~ h 2 + f-(2fc)!9n 2fc -! 6 ' { * } 

n=l u k=l v 1 n= Q 

where B 2k are the Bernoulli numbers, B 2k ~ {-l)^ 1 {2k)\ / ' {■n 2k 2 2k - 1 ). The 
leading order divergence comes from the integral, whereas the remaining 



Eigenvalue statistics in quantum ideal gases 



9 



terms give a power series in j3, 



zm ~ r /r 1/a + £ Cj p-^ . (29) 



a 

j 

A power series expansion of exp(— (3n a ) shows that only the powers 7 = 
(2k — 1)q carry non-zero weights. As noted before, these do not contribute 
to the density of states for the single particle but can be brought to live 
in the many particle system through combinations with sufficiently many 
powers of the first term. 

For rational a this does not work, since the derivatives required in the 
Euler-MacLaurin formula do not terminate, giving diverging coefficients for 
the powers of 0. Numerical results are consistent with the leading order 
behavior being given by the integral and the further terms in the desym- 
metrized version being a power series in E l l a . 

For later reference, I note the mean density of states for D particles that 
results from this asymptotic expansion, 

r («±n D 

p D (E) = V aj E(D - a)/a + . . . m 

The mean density of states thus decreases for D < a, is asymptotically 
constant for D = a and increases for D > a. 



3 Nearest neighbor spacings 
3.1 Numerical results 

The Maxwell-Boltzmann gas with its huge average degeneracies because of 
permutation symmetry gives rise to a rather singular level spacing distribu- 
tion, 

P(s) = (1 - g)S(s) + ±e~9 s (31) 

where the mean degeneracy g = D\. These massive degeneracies will ob- 
viously not survive perturbations. However, as Shnirelmans observations 
show, the low lying spectra of a weakly perturbed system may still have a 
c)-function contribution because of (almost) degeneracies between the per- 
turbed symmetry related tori (shnirelman 1993, Casati and Shepelyansky 
1995). 
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Within the symmetry reduced Bose-Einstein and Fermi-Dirac subspaces 
one does not expect this effect and generically finds a Poissonian spacing 
distribution (in accord with Berry and Tabor 1997a). However, for integer 
powers of a, all the energy levels will be integer and there can be some 
number theoretic degeneracies (Casati et al 1985). In the case of a = 2 it 
is known that the level spacing distribution collapses to a delta function. 
Thus although the density of states is constant, some large integers can be 
expressed in an increasing number of ways as sums of squares, and the gaps 
between them increase logarithmically. 

As to the case of cubes, the density of two cubes decreases, so that the 
unfolded levels (with mean spacing one) are not simple multiples of integers. 
Thus a continuous range of unfolded spacings can be achieved and the level 
spacing distribution is essentially Poissonian, except for a small overshoot 
at the origin (Fig. |2|a). 

For three cubes, the density is constant, and the spacing distribution seems 
to converge to 5-spikes with spacing l/p(E). The deviations noticeable in 
Fig. [2|b are due to finite E effects: the Bose-Einstein spacings have a wing 
to lower values because of the decreasing density, the Fermi-Dirac spacings 
have one to higher values because of the increasing density. However, for 
four and more cubes, a strong (5-function develops at the origin, eventually 
absorbing all the density (Fig. §c,d). 

3.2 Connection to Warings problem 

In the case of integer a, all eigenenergies are integer. Thus, if the density 
of states increases with energy, eventually one will reach a situation where 
the density exceeds one level per unit interval. This critical energy can be 
calculated to leading order from (^0|) to be 



For E > E c more and more levels have to fall onto the same integer, giving 
rise to a 5-function at the origin in the spacings distribution. This will 
happen for D > a + 1. The case D = a is marginal. For a = 2 it is 
known that the spacings between numbers that can be represented as sums of 
squares increases logarithmically, giving rise to a logarithmically increasing 
degeneracy: the spacing distribution converges to a delta function at the 




(32) 
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Table 1: Results on Warings numbers, taken from (Ribenboim 1991). The 
probabilistic lower estimate is often found to be too optimistic. The critical 



energies E c (computed using 32) in the last two columns increase rather 
rapidly. Some numerical consequences of this will be studied in section [O. 
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origin. For a = 3 and higher, the distribution seems to converge to a stick 
diagram. 

Since the density of states increases, one can ask whether all integers can 
in fact be represented as a sum of D integers raised to the power a. This 
is Warings problem (Ribenboim 1989). More precisely, define a number 
g{a) so that all integers can be represented if D > g(a). Since some small 
integers cause special problems, Waring considered another number G(a), 
such that if D > G(a) then all sufficiently large integers can be represented. 
Obviously g(a) > G(a) and G(a) > a + 1 because of the above density 
argument. Some results are collected in table 1. 

The considerations of the permutation symmetry give rise to a specializa- 
tion. Since both the Maxwell-Boltzmann case (which agrees with the sums 
Waring considered) and the Bose-Einstein statistics allow for repetition of 
integers, Warings numbers remain unchanged. However, the Fermi-Dirac 
statistics poses the additional constraint that all the integers have to be dif- 
ferent. Thus it is not possible to fill in l's if there is still a small gap to an 
integer. Therefore, the equivalent of g(a) makes no sense, some small inte- 
gers will always be missed. However, the density of states still increases for 
D > a + 1 without bound, so that the density of points is sufficient to reach 
all larger integers. Numerical tests suggest that at least for a = 3 and D = 4 
and 5 the density of points not represented among the lowest 10 6 integers 
decreases, but that is insufficient since the numbers involved rapidly grow 
large. It might be interesting to study the existence and values of G(a), 
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such that for D > G(a) all sufficiently large integers can be represented as 
sums of a's powers of D different integers. Obviously, G(a) > G(a). 

4 Form factor and pair correlations 

4.1 Numerical Results 

The level spacing distribution P(s) is a complicated mixture of n-point cor- 
relation functions, n = 2,3,... and thus not accessible to a complete semi- 
classical analysis. Some progress can be made for the two point correlation 
function 

C(e) = (p(E + e/2)p(E-e/2))/(pf (33) 

and the derived quantities spectral rigidity and number variance (Berry 
1985, Seligman and Verbaarschot 1987, Verbaarschot 1987, Bohigas 1991). 
The Fourier transform of the correlation function (|33l), the form factor, is 
the absolute value square of the Fourier transform of the spectrum. Since 
the latter can be related to periodic obits via the Berry- Tabor (1977b) semi- 
classical expansion or Poisson summation, a combination of classical results 
(Hannay and Ozorio de Almeida 1984) and quantum information (Berry 
1985) can be used to estimate the form factor. 

4.2 Berry- Tabor expansion and Poisson summation formula 

The density of states in the Bose-Einstein and Fermi-Dirac subspaces is to 
leading order given by that for the Maxwell-Boltzmann case with corrections 



for energy levels where two or more quantum numbers coincide (cf. 18). The 
Berry- Tabor (1977b) semiclassical expansion in terms of classical periodic 
orbits for the full density is thus a superposition of the ones for D or fewer 
particles with Maxwell-Boltzmann statistics. Technically, the semiclassical 
expansion of Berry and Tabor reduces to a Poisson summation on the EBK 
quantized eigenvalues, 



where 



( E -J2 n i) = ^g(mi,...,m D ) 

rti \ i=l / mi 

g(mi, • • • , m D ) = J d D xUg 1 Q(^ (x ~ £ ■■ 



(34) 



(35) 
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and 

r i x > o 

0(x) = < 1/2 x = (36) 
[ x < 

is the Heaviside step function. The Fourier transform of the product of step 
functions gives rise to another set of corrections to the leading order term, 
which can be combined with the ones due to symmetrization. 

Going through the algebra of expressing the delta function by its Fourier 
representation, doing the x-integrals and then the final /c-integral in sta- 
tionary phase, one ends up with an expression of the form (Seligman and 
Verbaarschot 1987) 

g(m u m 2 ,...,m D ) ~ {Um .y^-m^-i)) L (a-i-D)/(2a) . 

E -(2a-D-l)/(2a) gxp ^iL^^E 1 ^) (37) 

where 

L = Em W(-D. (38) 

Semiclassically, the L(mi, . . . , mu) are the periods of the classical motions. 
The lower order terms where one or two quantum numbers are equal are of 
similar form, but with coefficients multiplying the powers of m,. 

The main conclusion to be drawn from this is that the exponents a and 
a = a /{a. — 1) are conjugate to one another. In particular, a = cor- 
responds to a = Thus the periods for the eigenvalues with a = 3/2 
are sums of cubes. Since sums of four or more cubes show large degenera- 
cies, the periodic orbit spectrum for this system is degenerate which should 
influence the pair correlation function. 

4.3 Pair correlations 

By the Berry-Tabor (1977b) expansion the density of states can be written 

p(E) = J2 9m e^ El/a (39) 

m 

with the actions 

S m = L(m)( Q - 1 )/ Q (40) 

and the amplitudes given above. The Fourier transform in the scaling 
variable x = E l l a (alternativly one can expand around a reference en- 
ergy Eq — > oo and consider a small interval around Eq) then consists of 
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5- functions at actions S m . When taking the absolute value square, there are 
two contributions, one from every orbit with itself and one from different 
orbits. If there are no degeneracies in orbit actions, the diagonal part gives 
the constant form factor expected for Poissonian distributed levels (Berry 
1985). However, with degeneracies as in the case a = 3/2, there are addi- 
tional contributions from the cross terms and the form factor is higher by a 
factor g, the mean degeneracy of actions (Biswas et al 1991). 

The diagonal approximation can be good at best up to a period Tjj ~ 
Hp(E), for then the individual eigenvalues can be resolved and a quantum 
sum rule predicts the saturation to a constant, the Fourier transform of a 
5- function (Berry 1985). Since the degeneracy increases slowly with the pe- 
riod, its effect will only be noticeable if the density of states is sufficiently 
high so that for orbits of period near Th the degeneracies will become im- 
portant. For the lowest 100000 states studied here and in the studies of 
Biswas et al (1991), no effect of periodic orbit degeneracies on pair corre- 
lation functions was found. This may have to do with the large numbers 
involved before the mean density of states is so high that degeneracies are 
enforced. This problem awaits further study. 
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Figure 2: Level spacing distribution P(s) for a = 3 and various D. The 
dashed curves refer to the FD, the dotted curves to the BE subspaces. (a) 
D = 2, (h) D = 3, (c) D = 4 and (d) D = 5. The number of levels included 
in each diagram was about 96000 for each symmetry subspace. In (c) and 
(d) there is a strong delta function at the origin. The remaining features 
seen are not stationary and disappear as the number of levels is increased. 



